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Abstract. An existence and uniqueness result, up to fattening, for a 
class of crystalline mean curvature flows with natural mobility is proved. 
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1. Introduction 

In this paper we deal with the anisotropic mean curvature motion; that is, flows 
of sets 1 1 —>■ E(t) (formally) governed by the law 

(1.1) V(x,t) = -m(v E(t) )K E(t) (x), 

where V (x, t) stands for the (outer) normal velocity of the boundary dE(t) at x, <f> is 
a given norm on representing the surface tension , k E ^ is the anisotropic mean 
curvature of dE(t) associated with the anisotropy <f>, and to is a positive mobility 
which depends on the outer unit normal v E< P> to dE(t). Such an evolution law 
may be regarded as the gradient flow (with respect to a suitable formal Riemannian 
structure) of the anisotropic perimeter functional 

(1.2) P<t>{E) = [ </>(v E )dU N ~\ 

JdE 

the anisotropic curvature n E of dE being nothing but the first variation of (11.21) at 
E. When <j> is differentiable in 1^ \ {0}, then k e is given by 

(1.3) k e = div (V^)) . 

However, we are particularly interested in the case when is not differentiable, for 
instance the crystalline case, when the unit ball B$ := {</> < 1}, known as the Frank 
diagram , is a convex polytope. In the latter case, we will only consider the natural 
mobility given by m = <f. With this choice, (ED has the interesting property that 
the flow starting from a Wulff shape, that is, a level set of the polar </>° of </>, consists 
in a one-parameter family of shrinking Wulff shapes that extinguish in finite time. 
We recall that Wulff shapes are the only solutions to the isoperimetric problem 
associated with P ^ (see j26j). 
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The law (11.11) is used to describe several phenomena in Materials Science and 
Crystal Growth, see for instance :3Dj 32]. From the mathematical point of view, 
the geometric motion is well defined in a classical sense in the smooth case, that is, 
when (j> is at least C 3, “ (as well as the initial surface, and except at the origin) and 
“elliptic” (which means for instance that (j) 2 is strongly convex) [2|. Of course, the 
classical mean curvature flow falls within this class and corresponds to the choice 
(j) = Euclidean norm. In the smooth case, the main mathematical difficulties are 
related to the fact that singularities (like pinching) may form in finite time (see for 
instance [32]) in dimensions N > 3. Thus, the strong formulation of in]), which 
requires smoothness of the evolving sets, is well defined only for short times and 
one needs a weaker notion of solution that can handle the presence of singularities 
in order to define the flow for all positive times. When (j) is smooth, this task 
has been already accomplished and different approaches have been proposed in 
the literature, starting from the pioneering work by Brakke |15j . who suggested a 
weak formulation of the motion by mean curvature yielding deep regularity results 
but lacking uniqueness. These uniqueness issues have been subsequently overcome 
via the so-called level set approach [38] [MJ 121] . In particular, the case of cni) 
for m,4> of class C 2 is covered by [21] . The main idea is to represent the initial 
set as the zero sublevel set of a function uq and then to let all these level sets 
evolve according to the same geometric law (which makes sense thanks to the fact 
that the evolutions which we consider preserve inclusion). This procedure defines 
a time-dependent function u(x,t ) and transforms the geometric equation into a 
(degenerate) parabolic equation for u , which is shown to admit a unique viscosity 
solution with the prescribed initial datum uq. The evolution of the zero sublevel 
set of such a solution defines a generalized motion (see also [llj). which exists for 
all times and agrees with the classical one for short times, before the appearance of 
singularities (see [25]). Such a motion satisfies a comparison principle and is unique 
whenever the level sets of u have zero Lebesgue measure, i.e., whenever the so-called 
fattening phenomenon does not occurs. Fattening may in fact appear even for a 
smooth initial datum E° (see 0 ), but its occurrence is in some sense very “rare”: 
for instance, it is easy to understand that almost all the sublevels sets of the signed 
distance function from any given set E° will not generate any fattening. 

A third approach is represented by the minimizing movements scheme devised 
by Almgren, Taylor and Wang [2] and Luckhaus and Sturzenhecker [36]. It con¬ 
sists in constructing a sequence of discrete-in-time evolutions by iteratively solving 
suitable incremental minimum problems. Any limit of these evolutions as the time 
step vanishes defines a motion, which exists for all positive times (and is shown 
to be Holder-continuous in time for the L 1 norm). The connections between the 
generalized level set motion and Brakke solutions has been investigated in [34] . A 
simple proof of convergence of the Almgrem-Taylor-Wang (ATW) to the general¬ 
ized motion is shown in [20] . while a consistency result was already shown in [2]. 
See also [22] for a similar convergence proof in a more general setting (allowing for 
unbounded surfaces, as in the present paper), and [35] for new proofs and a gen¬ 
eralization to partitions. Roughly speaking, it turns out that whenever fattening 
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does not occur, the generalized level set motion coincides with the ATW flow and 
is also a solution in the sense of Brakke. 

Let us now consider a crystalline anisotropy. This case is more difficult, due to 
the lack of smoothness in the involved differential operators. Indeed, the crystalline 
normal S7(j)(y E ) is not uniquely defined for some directions and one needs to look 
at suitable selections of the (multivalued) subdifferential map, that is, vector fields 
z : dE —»• such that z(x) £ dcj)(v E (x)) for a.e. x. If there exists an admissible 

field z with tangential divergence div r z in L 2 {dE), then the crystalline curvature is 
given by the tangential divergence of z , where div T z has minimal L 2 - norm among 
all admissible fields (see m ni)- In particular, the crystalline curvature has a 
nonlocal character. 

Showing (even local-in-time) existence and uniqueness for crystalline mean cur¬ 
vature flows is somewhat harder and still largely open. Only in dimension 2, the 
problem has been settled by developing a crystalline version of the viscosity ap¬ 
proach for the level-set equation, see [28]. If the initial set is itself an appropriate 
planar crystal, the evolution equation boils down to a system of ODEs which has 
been studied in many former works, see in particular [3] HJ El Eli (Ml, while existence 
and uniqueness of strong solutions for initial “regular” (in an appropriate sense) sets 
was shown recently in HZ], One advantage of the level-set approach of [28] is the 
ability to address much more general equations where the speed depends on the 
crystalline curvature and the normal in an non-linear way. 

In dimensions N > 3, the only general available notion of global-in-time solution 
we are aware of is the minimizing movements motion provided by the ATW scheme; 
however, no general comparison results have been established so far. In fact, the 
higher-dimensional uniqueness results we know of deal with special classes of initial 
data (for instance convex initial data as in 16, [Hj or polyhedral sets as in [31]) 
or with very specific anisotropies (see [29] where a comparison principle valid in 
all dimensions has been established for the anisotropy </>(u) = |z/| + |i/jv|, with 
ujv := u-ejv and \v'\ the Euclidean norm of the orthogonal projection of v onto e^). 
However, Y. Giga has recently announced a very general existence and uniqueness 
result in the viscosity sense in dimension N = 3. 

In this paper we prove a global-in-time existence and uniqueness (up to possible 
fattening) result for the crystalline mean curvature flow valid in all dimensions, for 
arbitrary (possibly unbounded) initial sets, and for general crystalline anisotropies 
</>, but under the particular choice m = <j> in m- We do so by providing a suitable 
weak formulation of the problem and then by showing that such a notion yields a 
comparison principle. We then implement a variant of the ATW scheme to establish 
an existence result. 

Le us describe our approach in more details. It is based on ideas of [391 E|- In 
order to motivate our formulation, let us assume for a moment that <f> is smooth 
and that 1 H>• Eft) is a regular flow obeying (11.111 . Set d(-,t) := dist(-, E(t)), where 
dist denotes the distance induced by the polar norm cj>° (see (12.11) and (12.31) below). 
Then it is easy to see that the time partial derivative dtd of d on dE{t) equals 
—V/(f>(y E ( t d), with V denoting the outer normal velocity of the moving boundary. 
On the other hand, this quantity V/(j)(v E ^ t ') is nothing else as the speed of the 
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moving boundary along the Cahn-Hoffmann normal \/cj)(iz E ^), see [551 IT3] . Thus, 
ED may be rewritten as 

d t d = = div(V</>(Vd)) on dE(t) = d{d(-,t) = 0}. 

(Here and throughout the paper V stands for the spatial gradient.) On the other 
hand, if we look at a positive s-level set of d, the (weighted) normal velocity of 
x G {dist(-,£) = s} equals the normal velocity of its projection y on dEft), which 
is given by the anisotropic curvature {y ) of dEft ) at y. Since (as long as the 

surfaces are smooth) 

K {d ( .,t)=s } (>) = div(V(/>(Vd))(x,t) < K^ (t) (y), 

we deduce that 

(1.4) dtd > div(V</>(Vd)) in {d > 0} 

as long as E(-) is nonempty. In words, the positive level sets of the distance function 
shrink with a velocity which is higher than that given by the anisotropic curvature, 
and thus they may be regarded as super-flows or supersolutions of the geometric 
motion. Analogously, setting d c (-,t) := dist(-, E c (t)), where E c stands for the 
complement of E, we have 

(1.5) d t d c > div(V0(Vd c )) in {d c > 0} 

as long as E c (-) is nonempty. We may conclude that a smooth flow t i-a Eft) of sets 
solves ED if and only if (11.41) and ED are satisfied. 

As already remarked before, when <f> is crystalline V^(Vd) may not be defined 
and must be replaced in general by a suitable selection of the subdifferential map, 
that is, by a vector-field z 6 L°°({d > 0};®^) such that z(x) G dcj)(X7d(x)) for a.e. 
x , where df> denotes the subdifferential of (f>. Any such z will be called admissible 
for d. 

The above discussion motivates the following weak formulation of the crystalline 
flow: we will say that a one-parameter family t H > Eft) of closed sets, satisfying 
suitable continuity properties (see Definition 12.11 below) is a weak supersolution of 
ED with initial datum E° if E{ 0) C E° and there exists a vector-field z, admissible 
for d, such that ED hold in the sense of distributions, with V</>(Vd) replaced by 
2 . We will say instead that t Eft) is a weak-subsolution of (11.11) if E{ 0) 2 E° 
and t i—^ ( E(t)) c is weak supersolution. Finally, we will say that t i —> Eft) is a 
weak solution if it is both a weak sub- and a supersolution (with initial datum E°). 
Mostly for technical reasons, we will require in addition that the positive part of 
div 2 is bounded in {d > <5} for all S > 0. 

Let us notice that this formulation of the curvature flow in terms of the distance 
function has been already exploited for the standard mean curvature motion and its 
regular anisotropic variants. In fact, it is close in spirit to the distance formulation 
proposed and studied in (5H|, although it is somewhat stronger as it require the 
differential inequalities to hold in a distributional sense, rather than in the viscosity 
sense considered in |B9j . In this respect, our formulation is reminiscent of the 
approach developed in m- 
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We now describe the plan of the paper. In Section [51 after recalling some pre¬ 
liminaries definitions and introducing the main notation, we give the precise weak 
formulation of the sub- and supersolutions to the anisotropic mean curvature flow. 
In Section [3] we establish a comparison principle between sub- and supersolutions, 
which by standard arguments yields the uniqueness of the crystalline flow whenever 
fattening does not occur. We remark that the distributional formulation described 
above allows for a proof of the comparison, which is closer in spirit to the unique¬ 
ness proofs for standard parabolic equations. In particular, our argument is more 
elementary than the typical “viscosity” proof that is based on delicate regulariza¬ 
tion procedures and fine differentiability properties of semiconvex functions. In Sec¬ 
tion 2] we provide an existence results for the the weak formulation of the crystalline 
flow, which is based on the reformulation of the minimizing movements scheme of 
Almgren-Taylor-Wang / Luckhaus-Sturzenhecker introduced in [181 i 16]. Such a 
variant can be considered as a combination of the ideas of [2] and the threshold dy¬ 
namics algorithm studied in |23j . and has several advantages: for instance, it makes 
it easier to establish a comparison principle for the discrete-in-time evolutions and 
it works equally well for bounded and unbounded sets (as already exploited in [22]). 
In the main theorem of the section, we establish the convergence of the minimizing 
movements scheme to a weak solution, whenever no fattening occurs. 

We conclude this introduction by commenting on the restriction m = (f> in (O). 
Although such a mobility is rather natural (for instance it forces Wulff shapes to 
evolve in a self-similar way), it is not the most general case and different mobili¬ 
ties could be considered as physically interesting. However, at the moment, in the 
crystalline case we are able to provide the right convergence estimates for the mini¬ 
mizing movements scheme only under this assumption; the main technical reason is 
related to the fact that if dist is the distance induced by the polar norm (ff , then the 
crystalline curvatures of the positive level sets of dist(-, E ) are bounded above (this 
can be easily understood since in this case the sublevel sets of dist(-,.E) admit an 
inner tangent Wulff shape at all points of the boundary). Nevertheless, we remark 
that in the case of a smooth elliptic anisotropy, all our results and methods would 
work with any mobility m, thus showing that the viscosity solutions already studied 
in |2! ED Ell] satisfy in fact a stronger (distributional) formulation. The extension 
of our results to more general mobilities in the crystalline case will be the subject 
of future investigations. 
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2. A WEAK FORMULATION OF THE CRYSTALLINE MEAN CURVATURE FLOW 

In this section we introduce a suitable weak formulation of the crystalline mean 
curvature flow. Such a notion of solution resembles the formulation due to [39] . 
However, here we will not consider the viscosity setting of [39] and we will rather 
be concerned with distributional solutions (which appear for instance in 16j). 

2.1. Preliminaries. In this subsection we introduce the main objects and notation 
used throughout the paper. 

Let (j> denote a fixed norm on M. N , that is, a convex, even and 1-homogeneous 
real-valued function, which will play the role of the anisotropic interfacial energy 
density. In the terminology of crystal growth this is also called surface tension. 
Note that we do not assume any further regularity on (f> and in fact the main case 
of interest is when <j> is crystalline , that is, when the associated unit ball is a convex 
polytope. The interfacial energy is then given by 

P^E) :=sup{/ divC dx :(£Ct(R N ;R N ), </>°(C) < 1 

where we recall that the polar norm <j>° is defined as 

(2.1) <£°(0 := sup r)-£. 

<K?7)<1 

It can be checked that P^,{E) is finite if and only if A is a set of finite perimeter 
and, in this case, 

P<t>{E) = f ^ E )dn N ~\ 

J d*E 

where d*E denotes the so-called reduced boundary of E (see for instance [5j). More 
generally, given a function u £ BVt oc (M. N ) we may consider the anisotropic total 
variation maesure of u , which on the open (bounded if u BV(M. N )) subsets 
n C is defined as 

4b(Du)(fl ) := su P {/ udiv^dx : f £ (^(fijK^), </>°(C) < 1 

Because of the homogeneity of <j> it turns out that 4>{Du) coincides with the nonneg¬ 
ative Radon measure in R w given by Vudx + (f> \D s u\, where Vu stands for 

the absolutely continuous part of Du and | denotes the Radon-Nykodim deriv¬ 
ative of the singular part D s u of Du with respect to its (isotropic) total variation 
\D s u\, see [5]. 

Among the important properties of 4> and <p° let us mention the fact that d<j>( 0) = 
{£ : </>°(£) < 1} while d(j>°( 0) = {£ : </>(£) < 1}. Moreover, for rj ^ 0 

(2.2) d<j>{ri) = {£ : <j>°{£) < 1 and £ • ry = 4>{r])} = {£ : </>°(£) = 1 and £ • rj = 

(and the symmetric statement for (j>°). An easy consequence of the above charac¬ 
terization is that if r] £ d(f°{x) and x ^ 0, then x/cj> 0 {x ) £ dtfijj). 

The set 

IT(0,1) := {y : <t>°{y) < 1} 

is called the Wulff shape associated with </>. More generally, for x £ R w and R > 0, 
we will denote by 


W{x,R) := {y : <j>°(y~x) < R} 
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the Wulff shape of radius R and center x. In the Finsler metric framework associ¬ 
ated with (jf , Wulff shapes play the same role as standard balls do in the Euclidean 
setting. In particular, it is well-known that W (0, R) is the unique (up to transla¬ 
tions) solution of the anisotropic isoperimetric problem 

min {P^E): \E\ = | W(0, i?)|} , 

see for instance [25] , 

Given a set E C K^, we denote by dist(-, E) the distance from E induced by cj)° , 
that is, for any x £ 

(2.3) dist(x, E) := inf (f>°(x^y ) 

y£E 

if E / I and dist(x,0) := +oo. Moreover, we denote by d_E the signed distance 
from E induced by <f>°, i.e., 

dE{x) := dist(x, E) — dist(x, E c ) 

so that dist(x,E) = ds(x) + and dist(x, E c ) = ds(x)~ (here and throughout the 
paper we adopt the standard notation t + := t V 0 and t~ := (— t) + ). Note that 
</>(Vdu) = 1 a.e. in R^ \ dE. 

We finally recall the notion of Kuratowski convergence. We say that a sequence 
of closed sets E n in R m converges to a closed set E in the Kuratowki sense, and we 
write 

E n ^E, 

if the following conditions are satisfied: 

(i) if x n £ E n , any limit point of {x n } belongs to E\ 

(ii) any x £ E is the limit of a sequence {x n }, with x„ £ E n . 

One can easily see that E n —> E if and only if dist(-,£'„) —> dist(-,E) locally 
uniformly in R m (here one may consider the distance associated to any norm). 
In particular, by the Ascoli-Arzela Theorem, any sequence of closed sets admits a 
subsequence which converges in the Kuratowski sense. 

2.2. A weak formulation of the crystalline flow. In this subsection we intro¬ 
duce the weak formulation of the crystalline flow we will deal with. We refer the 
reader to the introduction for the motivation behind this definition. 

Definition 2.1. Let E° C R^ be a closed set. Let E be a closed set in M. N x [0, +oo) 
and for each t > 0 denote Eft) := {x £ R w : (x, t) £ E}. We say that E is a 
supersolution of the curvature flow cm with initial datum E° if 

(a) E{ 0) C E°; 

(b) for all t > 0 if Eft) = 0, then E(s) = 0 for all s > t; 

(c) E(s) Eft) as s /* t for all t > 0 (left-continuity); 

(d) setting d(x, t) := dist(x, Eft)) for (x, t) £ x (0, T*) \ E and 

T* := inf{t > 0 : E(s) = 0 for s>t}, 
then the inequality 


( 2 . 4 ) 


d t d> div^ 
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holds in the distributional sense in M. N x (0 ,T*) \ E for a suitable z G 
L°°(EL n x (0,T*)) such that z G dcj)fVd) a.e., divz is a Radon measure in 
R n x (0 ,T*)\E, and (divz)+ G L°°({(x,t) G R N x (0,T*) : d{x,t) > 5}) 
for every S > 0. 

We say that A, open set in x [0, +oo), is a subsolution with initial datum E° 
if A c is a supersolution with initial datum ( E°) c . 

Finally, we say that E, closed set in R w x [0,+oo), is a solution with initial 
datum E° if it is a supersolution and if E is a subsolution, both with initial datum 
E°. 

Remark 2.2. Notice that the initial condition for supersolutions may be rewritten 
as E° C A(0). In particular, if dE° = dE° and E is a solution according to the 
previous definition, then E( 0) = E°. 

Remark 2.3. If cf> is C 2 , then one can check that this definition is stronger than the 
definition in the viscosity sense (see in particular [39', j lT] ). 

We start by observing some useful continuity properties of the map d introduced 
in the previous definition. 

Lemma 2.4. Let E be a supersolution. Then, for eacht G [0,T*) ; df,s) converges 
locally uniformly in {x : d(x, t) > 0} as s \ t to for some function d r with d r > 
df,t) in {x : d(x,t) > 0}. 

Remark 2.5. Observe that by condition (c) in the definition (which is mostly tech¬ 
nical and forbids artificial constructions such as a supersolution which jumps to 
Eft) = 1^ at a given time i > 0), i i—X df,t) := df,Eft)) is left-continuous with 
respect to the local uniform convergence. 

Proof. By condition (d) ofDefinition l2.il the distributional derivative dtd is a Radon 
measure in R w x (0, T*)\E , so that d is locally a function with bounded variation 
in this (open) domain. In particular, for a.e. x G the map s i-A d(x , s) has a 
right limit d r (x, t) at each time t G [0, T*) such that d(x, t ) > 0. Since the functions 
d(-,s ) are also equi-Lipschitz in space as s varies, we may conclude that the right 
limit is in fact locally uniform in {x : d{x,t ) > 0}. 

Since E is closed, for every t G [0, T*] we clearly have that all Kuratowski 
cluster points of E(s) as s —> t are contained in Eft), equivalently, d(x,t) < 
liminfg-K dfx, s). Thus, d r > df,t) in {x : d(x,t) >0}. □ 


3. Comparison results 


In this section we prove the main comparison principle between sub- and super¬ 
solutions (see Theorem 13.31) . In Lemma EOl below, we establish a first (suboptimal) 
comparison result between a supersolution and a suitable anisotropic total variation 
flow (see HEZl). To this aim, we give an explicit solution to the anisotropic total 
variation flow with initial datum cf>°. 


Lemma 3.1. The pair (/, <f) defined by 


(3.1) 


f{x,t) 


r(t) + t^jjf if (j)° (x) < rft) := y/(N + l)t, 

( a; ) + t IpTT) otherwise 
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and 

if <t>°{x) < r(t), 
if > r(t), 

solve the following Cauchy problem for the cf-total variation flow in R w : 

) d t f = div£ a.e. in R N x (0, +oo), 

0 £ <90(V/) a.e. in x (0,+oo), 

/(■) 0) = (f>°. 

Moreover, given A > 1, the pair (f\,(\) given by 

fx(x, t ) := A/O, t/X) Cx(x, t) := C(x, t/X) 

for (x,t) £ R w x (0, +oo) solves (13.31) . with the initial datum <f>° replaced by A0°. 

Proof. Recalling that / £ <90(V/) is equivalent to 0°(£) < 1, / • V/ = 0(V/) (see 
m), the proof follows by direct verification. The details are left to the reader. □ 

Next lemma provides a first comparison estimate, which is far from being sharp. 
However, the optimal estimate can be established a posteriori as a consequence of 
our main comparison theorem (see Theorem 13.31 below 1. 

Lemma 3.2. Let E be a supersolution and d := dist(-, E(-)) the associated one 
parameter family of distance functions. Assume that for some (. x , t) £ BO x [0, +oo) 
we have d{x , t) > R > 0. Then, there exists a constant Xi v > 0 such that d{x, t+s) > 
R-XNy/s for all s £ [0, R 2 /(16y^)]. 

Proof. Observe first that thanks to Lemma \‘2. 41 since d{-,t) > R/ 4 in {x : (f>°(x — 
x) < 3i?/4} = W{x, 3i?/4), there exists a (unknown) time t* such that d(-, t + s) > 
a > 0 in lT(x,3i?/4) for all s £ [0,f*] for some positive a. We will compare d with 
the solution S of the ^-total variation flow starting from 

(5(-,0) := R— ~0°(■ - x). 

More precisely, setting S(x, s ) := R— fi/z(x — x, s), where fi/ 3 {x, t) := 4/3 f(x, 3t/4) 
and / is given by m, by Lemma T3. II S satisfies 

dtS = div^ in R N x (0, +oo), 

£ £ <90(V<5) a.e. in R^ x (0,+oo), 

where t;{x,t) = — 0(x,3t/4), with / defined by (13.21) . Note that 6 is negative outside 
W(x,3R/4) for all positive times. 

Let T(s) be a smooth, convex, nonnegative function, which vanishes only for 
s < 0, and consider the function w(x,s ) := (<5(cc, s) — d(x,t + s)). Without loss 

of generality, we assume to simplify the notation that t = 0. By construction, 
w{x, 0) = 0 in W(x, 3i?/4) and w(-,s) = 0 on dW(x , 3i?/4) for 0 < s < t*. 

Since 0(Vd) < 1 a.e. and dtd is a measure wherever it is positive, it follows that 
d is a function in BVi oc (W(x , 3R/ 4) x (0 ,t*)) and its distributional time derivative 
has the form 

d t d = ^^[d(-, t + 0) — d(-, t — 0)]dx + dfd 

t£j 



(3.2) 


C(x,t) := <j ^ 

4 >° (x) 



10 ANTONIN CHAMBOLLE, MASSIMILIANO MORINI, AND MARCELLO PONSIGLIONE 


where J is the (countable) set of times where d jumps and dfd is the diffuse 
(Cantor+absolutely continuous) part of the derivative. It turns out that ('Lemma |2.4l) 
d(-,t + 0) — d(-,t — 0) >0 for each t e J. Moreover, since the positive part of di vz 
is absolutely continuous with respect to the Lebesque measure, (El entails 

dfd> divz. 

Using the chain rule for BV functions, see 0 ), one has 

d t w = ^[\H(5(-,t) - d(-,t + 0 )) - \H(<S(-,t) - d(-,t - 0 ))]dac 
teJ 

+ *'(6 - d){d t S - dfd) < - d)(div£ - divz). 

Hence, for a.e. t < t* , using the fact that (f> and T are convex, (<5 — d) vanishes 
on dW(x, 3R/4) and recalling (13.41) . we have 

dt f wdx < j 4''(<5 — d)(div£ — divz) 

JW(x,3R/4) JW{x,3R/4) 

= - [ (£-z)-(V<5-Vd)^"(d-d)< 0. 

Jw(x,3R/4) 

It follows that w = 'if (S — d) = 0, that is, d > S a.e. at all times less than t*. More 
precisely, for 0 < s < t* we have 

4 N I I 

(3.5) d(x, t + s)> R - f 4/3 {x - x, s) = R - N + 1 =■ R~ XnVs. 

It follows from (13.51) that d(x,t + s) > 3R/4 and, in turn, d(-,t + s) > 0 on 
dW{x, 3R/4) for all s < min{t*, I? 2 /(16y^)}. But then we can restart the ar¬ 
gument above to find that (13.51) remains valid for slightly larger times. Thus, we 
may conclude that El holds at least for all 0 < s < R 2 /(16 xn)- This concludes 
the proof of the lemma. 

□ 

Now we can state the main result of this section, which is a comparison result 
between sub- and supersolutions. 

Theorem 3.3. Let E be a supersolution with initial datum E° and F be a subso¬ 
lution with initial datum F°. Assume that dist (E°,F° C ) =: A > 0. Then for each 
t> 0, dist (E(t), F c (t)) > A. 

Proof. Let TJ. and Tp be the maximal existence time for E and F. For all t > 
min {Tp,Tp} we have that either E or F c is empty. In this case, clearly the con¬ 
clusion holds true. 

Now, consider the case t < min {Tp,Tp} (and assume without lost of generality 
that Tp,Tp > 0). Let us fix 0 < r]i < < rf{ < rj'f < rf 2 < 772 < A. We 

will show the conclusion of the theorem for a time interval ( 0 ,i*) for a suitable t* 
depending only on 771 , rj[, 77 ", rj'f, vf ^, 772 , and ultimately only on A. It is clear 
then that reiterating the argument yields the conclusion of the theorem for all 
times. We recall that dp(x,t) := dist (x,E(t)) — dist (a :,E c (t)) and g?f is defined 
analogously. We denote by Ze and zp the fields appearing in the definition of 
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super- and subsolutions (see Definition ^. Ill , corresponding to E and F, respectively. 
Define 


S := {x £ : rj\ < <1e{x , 0) < 772 } 

and note that by Lemma EJT21 there exists t* > 0 depending only on 771 , A 


that 

(3.6) 


d E {x, t) > d E (x, 0 ) - xjv\/t 
d F {x, t) < d F {x, 0 ) + XN\/t 


for all x £ S and t £ (0,f*). 


772 such 


We now set 


d E '■= d E V (?7i + XnVI) , 
d F := (d F + A) A (r]' 2 - XnVI) . 
Clearly, by our assumptions d E (-,0) > d F (-, 0). We claim that 
(3.7) d F > d F on dS x (0 ,t*). 


Here and in the rest of the proof we may assume without loss of generality that t* is 
as small as needed (but still depending only on A). To this aim, write dS = TiUl^, 
where Id := (-, 0) = 771 } and T 2 := {g(e(-,0) = 772 }. Since d F (-, 0) + A < 

d E {-, 0 ) = 771 on Li, we deduce 

d F < d F + A < 771 + xnVI < rj[ < d E 


on Ti x (0, £*). Similarly one can show that the inequality d E > d F holds on 

r 2 x (0,**)■ 

Again by (13.61) we have 

(3.8) dE>y>0 in{d B (-, 0 )> 77 "}x (0,t*) 
and, observing that d F (-, 0) < 772 — A in {d£?(-, 0) < 772 }, 

(3.9) d F <— y—< 0 in {^(-,0) < 773} x (0,f). 

In particular 

E(t) CC F(t) for t £ (0, t*). 


We now claim that, setting 


S" := {x £ : r{\ < d E {x, 0 ) < 772 }, 

we have 


(3.10) d E = d E and dp = dp + A in S" x (0, t*). 

Indeed by (13.61) we have 

d E (x, t ) > 77 " - XnVI > 771 + xnVI for (x, t ) £ S" x ( 0 , t*) 


and thus d E = d E in S” x (0 ,t*). The proof of the second identity in (13.101) is 
analogous. 

Now we will use quite standard parabolic maximum principles, like in the proof 
of Lemma 13.21 Notice that 


d t d E = ^2[d E {-,t + 0 ) - d E (-,t — 0 )]dx + dfd E , 

tpj 
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where J is the (countable) set of times where d E possibly jumps and dfd E is the 
diffuse part of the distributional derivative. Using for instance the chain rule proved 
in j4], in S' x (0, t*) we have that 


dfd E = 


- ) 2Vt a - e - in : + XnVI > d E (x )} . 


dtd E \d?d E \-a.e. in {(x,t) : p[ + XnVI < d E {x)} . 


An analogous formula holds for d?d F . Recalling that (div 2 e) + and (di vz F ) belong 
to L°°(S x (0, £*)) it follows that (possibly modifying t*) 

(3.11) dfd E > di vz E and d?d F < divzp 

in the sense of measures in S x (0,£*). Note also that a.e. in S x (0,£*) 

(3.12) z E £ d(j){yd E ) and zf £ d(j)(Vd F ). 

Fix p > N and set 'F(s) := (s + ) p and w := \H( d F — d E )■ By (13.711 we have 

(3.13) w; = 0 on dS x (0, £*). 

Using as before the chain rule for BV functions, recalling (13.111) and the fact that 
the jump parts of dtd E and dtd E are nonnegative and nonpositive, respectively, we 
have 

(3.14) d t w < d>'(d F ~ d E )(dfd F — dfd E ) < ’F , ( d F — d E )(dxvz F — divz E ) 


in S' x (0, £*). Choose a cut-off function rj £ C%°( K”) such that 0 < rj < 1 and rj = 1 
on B\. For every e > 0 we set p E (x) := p(ex). Using (13.131) and (13.141) . we have 

d t [ wp p dx < f p p d!'(d F — d E )(d\vz F — d\vz E ) 

Js Js 

= — [ p E ^"(d F — d E )(z F — z E ) ■ (Vd F — Vds) dx+ 

Js 

p [ ry p_1 d>'(d F - d E )Vp s ■ (z F - Z E ) dx 

Js 

<P f Ve _1 ^'(d F - d E )Vp e ■ (■ Z F ~ Z E ) dx , 

Js 

where we have also used the inequality (z F — z E ) ■ ( Vd F — Vd E ) > 0, which follows 
from (13.121) and the convexity of (j). By Holder Inequality and using the explicit 
expression of and 'F / , we get 

r _ //■ x 1 - 1 

d t 


^wrf e dx<Cp HVrfclU^) \^J^wp p dx 

for some constant C > 0 depending only on the L°°-norms of z E and z F . Since 
u> = 0att = 0,a simple ODE argument then yields 

[ wrf e dx < (C'p||V77 e || L p (R jv ) £) p 

Js 

for all t £ (0,£*). Observing that ||V%|| p p(RiV) = e p_JV ||Vr/||^ p(RAr) -A 0 and p e Z' 1 
as £ —> 0+, we conclude that w = 0, and in turn d E > d F in Sx (0, £*). In particular, 
by claim (13.101) . we have shown that d E > d F + A in S" x (0, £*). We finally claim 
that dist(E(t), F c {t)) > A for t £ (0, t*). To see this, fix e > 0, and let let x £ dE(t) 
and y £ dF(t) be such that cj)°(x — y) < dist(U(£), F c (t)) + s. Note that by (13.81) and 
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(ED we have d E (i r,0) < rff and d E (y, 0) > rff. Thus there exists z € S" D [x,y\, 
where [a :,y] denotes the segment joining x and y. Since (-,£) > dp{-,t) + A in 
S", we have 

(3.15) dist (E(t), F c (t)) > (f°{x — y) — e = (f°(x — z) + <j>°(z — y) — e > 

— d,F(z,t) + d E (z,t ) — £ > A — £. 

The claim follows by the arbitrariness of e, and this concludes the proof of the 
theorem. □ 


4. Existence via minimizing movements 


In this section we prove an existence result for the crystalline curvature flow, 
according to Definition ^. II Such a solution is obtained via a variant of the Almgren- 
Taylor-Wang minimizing movements scheme ([2]) introduced in mm- 


4.1. Minimizing movements. Let E° C be closed. Fix a time-step h > 0 
and set E® = E°. We then inductively define E k+1 (for all k G N) according to the 
following procedure: If E k ^ 0, R N , then let {u k+1 ,z k+1 ) : R^ ->lx R N satisfy 


—h divz^ +1 + u k h +1 = d E k , 
z k h +1 G dcj){\7u k h +1 ) a.eAn R", 


and set E^ +1 := {x : u^ +1 < 0}. If either E k = 0 or E k = R^, then set E^ +1 := E k . 
We denote by T£ the first discrete time hk such that E k = 0, if such a time exists; 
otherwise we set = +oo. 

In proposition 14.11 below we will show that this construction is well defined, 
since problem ( 14 . 11 ) admits a unique solution u^ +1 that is Lipschitz continuous. In 
particular, E k+1 is a closed set for all k. 

Before stating the main facts about the differential problem ( 14 . 11 ) , we recall that 
given z G L°°{R N ;R N ) with divz G Lf oc (R N ) and w G BVi oc (R N )nLf oc (R N ), z-Dw 
denotes the Radon measure associated with the linear functional 


Lip := 

see [lUj . 



w pdivz dx 



vj z ■ Vpdx 


for all p G C' c °°( R n ), 


Proposition 4.1. Let g G Lf oc (R N ). There exists a field z G L°°(R N ; W(0,1)) and 
a unique function u G BVi oc (R N ) D Lf oc (R N ) such that the pair (u, z) satisfies 

{ —hdivz + u = g inT>'(R N ), 

4>°(z) <1 a.e. in R N , 

z ■ Du = (j>{Du) in the sense of measures. 

Moreover, for any R > 0 and v G BV(B R ) with Supp (u — v) <s Bn, 

<j>(Du)(B R ) + ^ J b ( u ~ d) 2 dx < </>(Dv){B R ) + ^ ( v-g) 2 dx , 

and for every s£l the set E s := {x G M N : u(x) < s} solves the minimization 
problem 

min P^F; B R ) + ^- I (g{x) - s ) dx. 
t /xpjs^Br ri Jfc\Br 
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If gi < g 2 and if u\, U 2 are the corresponding solutions to (14.21) (with g replaced 
by g 1 and g 2 , respectively), then u\ < 112 - 

Finally if in addition g is Lipschitz with </>(V g) < 1, then the unique solution u 
of (|4.2I) is also Lipschitz and satisfies (f>(Vu ) < 1 a.e. in R w . As a consequence, 
m is equivalent to 


(4.3) 


—h divz + u = g, 
z £ d(j){Wu) a.e. in 1^ 


Proof. See [TH Theorem 2], [TJ Theorem 3.3]. 


□ 


Remark 4.2 (Consistency with the ATW scheme). When dE° is bounded, the min¬ 
imality property of the level sets stated above shows, in particular, that the sets E£ 
are constructed according to the Almgren-Taylor-Wang scheme [2] . 


Since by the previous proposition (f>(Vu^ +1 ) < 1 a.e. in R N , one deduces, in 
particular, that 


(4.4) 


u 


k +1 


< d E k +1 in {x : dist(a;, £'^ +1 ) > 0} 


,fc+i 


u. h > d E k+i in {x : dist (x,E^ + ) < 0} . 

We are now in a position to define the time discrete evolutions. Precisely, we set 
E h := {( x,t) : x £ E l * /h] }, 


(4.5) 


E h (t) := E l £ /h] = {x : (x,t) £ E h }, 
d h (x,t) := d Eh ( t )(x), 

Uh(x,t ) := u^ h \x), 

z h (x,t) := z l * /h] (x), 


where [•] stands for the integer part of its argument. 


Remark 4.3 (Discrete comparison principle). The last part of Proposition l4Tl clearlv 
implies that the scheme is monotone, that is, the discrete evolutions satisfy the 
comparison principle. More precisely, if E° C F° are closed sets and if we denote 
by Eh and Fh the discrete evolutions with initial datum E° and F°, respectively, 
then E h C F h . 


4.2. Comparison with the Wulff shape. In this subsection, we exploit Re¬ 
mark [O] to compare the discrete evolutions (14.51) with the minimizing movements 
of the Wulff shape and derive an estimate, which will be useful in the convergence 
analysis. The evolution starting from a Wulff shape W(Q,R) is explicitly known. 
Indeed, from [16] Appendix B, Eq. (39)], the solution of (14.21) . with g replaced by 
dw(o,R) = (f>° - R, is given by (j)° h - R, where 


(4.6) 


<t>°h( x ) 


if - V h ( N + l), 

0 ° ( x ) + else - 


It follows that if E° = W(f),R), one has E h (t) = W(0,rf((t)) for a function that 
satisfies 


R+ ySR 2 - Ah(N - 1) 







EXISTENCE AND UNIQUENESS FOR A CRYSTALLINE MEAN CURVATURE FLOW 15 


if h < i? 2 /(4(7V + 1)). In particular, 

r%{h) > s/R 2 - Ah(N - 1) 

for the same h’s. By iteration, we have r^(t) > \JR 2 — 4 t{N — 1) > for 0 <t< 
R 2 /(8(N — 1)) and h < R 2 /(8(N + 1)). Since rff(t) = R for t £ [0, h), we infer 

(4.7) r*{t) > ^/R 2 -4t(N-l) 
for 0 < t < R 2 /(8(7V + 1)) and for all h. 

Now we return to the motion from an arbitrary set E°. If for some (x,t) £ 
M. N x [0,T£) we have dh{x,t ) > R, then W(x,R) fl Eh(t) = 0. Hence, by the 
comparison principle stated in Remark 14.31 and by (14.71) we have 

dh(x,s) > \JR 2 — 4(N — 1 )(s — t + h) 

for t < s and s + h — t < R 2 /(8 (N + 1)). 

By letting R /• dh(x, t ) we obtain 

(4.8) d h (x,s) > yJ d 2 h (x, t) - 4 (N - 1 )(s-t + h) 
for t < s and s + h — t < d 2 (x, t)/(8(N + 1)). 

4.3. Convergence of the scheme. Up to a subsequence we have 

E hl ^E and {E hl )° A A c 

for a suitable closed sets E and a suitable open set A C E. Define E(t) and A[t) 
as in (14.51) . 

Observe that if E(t) = 0 for some t > 0, then (14.81) implies that E(s) = 0 for all 
s > t so that we can define, as in Definition 12.11 the extinction time T* of E, and 
similarly the extinction time T'* of A c . Notice that at least one between T* and 
T'* is +oo. Possibly extracting a further subsequence, we have the following result: 

Proposition 4.4. There exists a countable setAf C (0,+oo) such that dh,(-,t) + 
dist(-, E(t)) and dh t (■, t)~ dist(-, A c ) locally uniformly for all t £ (0, +oo) \ Af. 

Moreover, for every x £ R w the functions dist(a;, £'(•)) and dist(a;,H c ) are left 
continuous and are right lower semicontinuous. Equivalently, the functions E(-) 
and A c are left continuous and are right upper semicontinuous with respect to the 
Kuratowski convergence. Finally, E{ 0) = E° and A(0) = E°. 

Proof. By the Ascoli-Arzela Theorem and a standard diagonal argument, we may 
extract a further (not relabeled) subsequence such that dh,(-,t) —> d(-,t) locally 
uniformly for all t. £ Q D (0,+oo), where d(-,t) is either a Lipschitz function or 
infinite everywhere. In the latter case, either d(-,t) = +oo or d(-,t) = — oo. 

We observe that for all t £ (0 ,T*) D Q we have d(-,t) < +oo. To see this we 
argue by contradiction assuming that for every x £ M. N and for every M > 0 we 
have dh t (x, t) > M for all l large enough. We may now apply (14.81) to deduce that 
there exists a right interval (f, t') independent of l such that dh t (x, s) > jf for l 
large enough and for all s £ that is, dhi(-,s) -4 +oo for all s £ (t,t'). This 

in turn would imply E(s) = 0 for all s £ (t, t’), which is impossible since t < T*. A 
similar argument shows that for all t £ (0 ,T'*) we have d(-,t) > —oo. 
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Let x £ R N and t > 0 be such that limsup; dh, (x, t) =: R > 0. Then, 
given 0 < R' < R, we have dh,(x,t) > R' for infinitely many l. By (14.81) . for 
t < s < t + CR' 2 we deduce limsup, dh, (x, s ) > R' 2 — 4 (AT — l)(s — t). In par¬ 

ticular if s £ Q fl (t,t + CR 2 ), then it follows d(x,s) > \/ R 2 — 4 (AT — l)(s — t). 
If, in addition, also t £ Q, then we have d(x, s) > \Jd(x,t) 2 — 4(7V — l)(s — t). 
Now let x £ and t > 0, and assume R := limsup se Q d(x, t) > 0. Con¬ 
sider a sequence of rational numbers Sk \ t such that lim*, d(x, Sk) = R. For 
s > t rational and close enough to t, if k is sufficiently large, then Sk < s and 
d(x, s ) > \Jd(x, Sk) 2 — 4 (N — l)(s — Sk)- Sending k to infinity it follows d(x, s) > 

yjK * 1 — 4 (N — l)(s — t) so that liminf sg Q j d(x, t) > R. Hence d(x, •) has a right 
limit at t. The same conclusion holds if lim inf seQ, a \t d(x, t) < 0, with the same 
proof. We deduce that the d admits a right limit (locally uniformly in space) at 
any t > 0. A similar argument shows that d also admits a left limit at any t > 0. 
Moreover, arguing similarly and using (14.81) again, we can show that 

d(x, t -(- 0) ± > lim sup dh, (x, s) 1 * 1 

Z—>-oo,s—>■£ 

> lim inf dh^XjS)^ > d(x,t — 0) ± . 

I —^oo, s—yt 

Let A f be the set of all times t such that the left and right limits of d differ at 
(x, t), for some x £ R w (we also assume 0 £ J\f). Notice that A f is countable, since 
it can be written as the union over k £ N and x £ Q N of the times such that the 
gap between the right and left limit of d(x, ■) is larger than 1 jk (which for fixed k 
and x cannot have cluster points). We denote by d(x,t) the common value of the 
right and left limits of d(x, ■) at t AT. 

By (|4.9D we immediately have that lim/^oo dh t (•, t) = d(-,t) for all t ^ AT. We now 
show that for t qL AT, we have d(-,t) + = dist(-, E(t)). This is equivalent to showing 
that E(t) coincides with the Kuratowski limit K of Eh, (t), since d(-,t) + = dist(-, K). 
Clearly, K C E(t). Conversely, if x ^ AT, then d(x,t) + =: R > 0. Since d is 
continuous at t , we may find e so small that lim;-^ dh, ( x , t — e)> d(x , t — e) > R/2 
and in turn, by (|4.8jl . W(x, R/ 4) x [t — e, t + e] D Eh, = 0 for l large enough. Thus 
x ^ E(t), showing that E(t) = K and d(x,t) + = dist (x,E(t)). A similar argument 
yields that d(x,t)~ = dist (a;, A c ). 

Always by (14.81) . one can easily prove that E( 0) C E°. Since Eh,(0) = E° for all 

l , we infer the equality E( 0) = E°. Symmetrically, one can show that A(0) = E°. 

Finally, we prove the continuity properties of E(t). The right upper semiconti¬ 
nuity with respect to the Kuratowski convergence is a consequence of the fact that 
E is closed. Let us prove now the left continuity. To this aim, denote by K the 
Kuratowski limit of E(s) as s /* t. Clearly K C E(t). Let now x ^ K. Then 
lim s /i t dist (a;, E(s)) = dist(a;,AT) =: R > 0. Arguing exactly as before we may 
choose e so small that liminf/dist(a;, Eh, (t — e)) > dist(a;,i?(t — e)) > R/2 and 
W(x, R/ 4) x [t — e,t + e] D Eh, = 0 for all l large enough, so that x £ E(t). Hence 
K = E(t). This establishes the Kuratowski left-continuity of E(-) and concludes 
the proof of the proposition. □ 

Theorem 4.5. E is a supersolution in the sense of DeRnition \ 2. 1\ with initial datum 
E°, while A is a subsolution with initial datum E°. 
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Proof. Points (a), (b) and (c) of Definition 12.11 follow from Proposition 14.41 It re¬ 
mains to show (d). Possibly extracting a further subsequence and setting Zf H (•, t) := 
0 for t > T ^ if T)* <T* 1 we may assume that Z} M converges weakly-* in L°°(R N x 
(0, T*)) to some vector-field z satisfying (p°(z) < 1 almost everywhere. Recall that 
by (14.411 we have u k+1 < d E k+ 1 , whenever d E k +1 > 0. In turn, it follows from (14.111 
that 

d E k+1 d E k 

(4.10) divz£ +1 < —-— -- a.e. on {d E k +1 > 0}. 

h h 

Consider a nonnegative test function y £ (^((“((IR^ x (0, T*)) \ E). If Z is large 
enough, then the distance of the support of rj from Eh t is bounded away from zero. 
In particular, dh t is finite and positive on Supp rj. We deduce from (14.101) that 


J j rj(x, t) ^ dh ^ x ' t + — d hl {x,t) _ diy2 ,^ ^ t + f ll ^j dtdcc 

^ x ' 1 — —dh t (x, t) - Zh, (x, t + hi) ■ Vr](x, t)^j dtdx > 0. 

Passing to the limit Z —>• oo we obtain (12.41) . 

Next, we establish an upper bound for divz^, away from Eh t ■ To this aim observe 
that 

d E k = min f (• - y) 
h 

so that, by (14.11) and the comparison principle stated at the end of Proposition l4.ll 



u 


fc+1 

h 


< min - y) 

yeK 


where (/% is given in (14.61) . Thus, if d E k{x) > R > 0, then 


,k +1 


(x) < min <j)°{x — y) + h 


N - 1 
R 


— d E k (x) -)- h 


N - 1 
R 


provided h < R 2 /{N + 1). As a consequence of (14.11) . we obtain 

(4.11) dwz k h +1 < a.e. in {x : d E k(x) > R}. 

it h 

It is then easy to deduce from the convergence properties of E'h, and dh t that 


divz < N in {(x,t) £ x (0,T*) : d(x,t) > R} 

R 

in the sense of distributions. It follows that di vz is a Radon measure in M. N x 
(0 ,T*)\E, and (divz)+ £ L°°({(x,t) £ R N x (0,T*) : d(x,t) > <5}) for every (5 > 0. 

We now provide a lower (Zi-dependent) bound for divzh, ■ To this aim, note that 
if d E k (x) =: R> 0, then d E k > R — <p°(- — x). Thus, by comparison as before, 


u k h +1 {x) > R 


cfKO ) = R-Vh 


2N 

VnTT' 


In turn, by (14.11) . we deduce 


di vz 


fc+i 

h 


1 2 N 

~ ~VhVN Tl 


a.e. in {x : d E k{x) > 0}. 
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Combining the above inequality with (14.1111 and using (14.11) again, we deduce that 
for all t £ (0, T*) \ J\f (where recall that J\f is introduced in Proposition 14.41) and 
any 6 > 0 

2 N 

\\ u h, {')t) d hl {-,t hi) \\L°°({x:dh l (x,t— h,)>6}) — » ^ ^ ’ 

provided that l is large enough. In particular, recalling the convergence properties 
of and dh, (see also (14.9D ~). we deduce that 

(4.12) u h , -> d a.e. in R N x (0, T*) \ E, 

with the sequence {uh,} locally (in space and time) uniformly bounded. 

Consider now, as before, a nonnegative test function r] £ C£°((R N x (0 ,T*))\E). 
Then, recalling (14.121) . we have by lower semicontinuity 

J j (f>(Vd)r) dxdt < lira inf J j (f>(Vuh, )r] dxdt = liminf J J \zh, • dxdt. 

On the other hand, 

(zh, ■ Vuh t )r)dxdt = / / (zh, ■ \7d)r]dxdt + / / Zh, ■ V(u? li — d)r]dxdt, 

with 

J J {zh t • Vd)r]dxdt ? J J (z ■ Vd)ry dxdt. 

Hence, we obtain 

(4.13) J J (j){S7d)r] dxdt < J J (z ■ Vd)r/ dxdt, 

provided we show that 

(4.14) lim J J Zh, -V(uh, — d)r]dxdt = 0. 

For each t, set 

mi(t) := min (uh,{x,t) — d[x,t )), MAt) := max (uh,(x,t) — d(x,t)). 

xGSupp ? 7 (-,t) aiGSupp ri(-,t) 

Recall that these quantities are uniformly bounded and converge to 0 at all t (jL AT. 
Then, we can write 

(4.15) J J z h , ■ V (uh, - d)r] dxdt = J J z h , ■ V (u h , - d - mi)r) dxdt 

= - J J (Uh, —d — mi){z h , ■ Vry + rjdivzh,) dxdt. 

For l large enough, since the support of 77 is at positive distance from E there exists 
<5 > 0 such that dh, > S everywhere on this support, so that div,^, < (N —1)/S. It 
follows that 


(uh, — d — mi)^dWzh, dxdt > — 


TV — 1 


(uh, — d — mi)r] dxdt 0, 


thanks also to (14.121) . Recalling (14.151) . we can conclude that 

liminf J J Zh, • V (uh, — d)rjdxdt > 0. 

In the same way, writing now 

J J z h , • V (u h , - d)r] dxdt = J J Zh, ■ V (uh, - d - M t )r] dxdt 
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and using u — d— Mi <0 a.e. on Supp 77 , one can show that 

limsup J j z h, ■ ^(uhi ~ d)ridxdt < 0 

so that (14.141) follows. In turn, (14.131) holds, that is, <j){S/d) < z ■ Vd a.e. in 
R w x (0, T*) \ E. On the other hand, recalling that cjf{z) < 1 a.e. in 1x (0,T*), 
we have 

z • Vd < (j){S7d) 

a.e. in x (0,T*). We conclude that cj){S7d ) = z ■ Vd and, in turn, z G <9</>(Vd) 
a.e. in R^ x (0 ,T*)\E. This concludes the proof that if is a supersolution. The 
proof that A is a subsolution is identical. □ 

Corollary 4.6. Let u° be a bounded, uniformly continuous in R w . Then for all 
s £ R but a countable number, the minimizing movement scheme starting from 
E° s = {m° < s} converges to the unique solution of the curvature flow in the sense 
of Defi,nition \2.R with initial datum E®. 

Proof. The arguments are standard and rely on the comparison theorem 13.31 The 
bad (countable set) is the set of levels for which “fattening” occurs, that is, |if\A| > 
0. Observe that from Theorem 13.31 one easily shows the existence of a unique level- 
set solution u(x,t) starting from u°, which shares the same spatial modulus of 
continuity and is also uniformly continuous in time (see for instance [19] Subsec¬ 
tion 6.3]). □ 


5. Conclusion and perspectives 

In this note we have shown the existence and uniqueness of a mean curvature 
flow (namely, the “natural” flow by mean curvature along the Cahn-Hoffmann vector 
field) with a technique which does not require any type of regularity on the surface 
tension, and thus have provided the first sound definition of a crystalline curvature 
flow in any dimension. It does not require that the initial surface is bounded and 
applies, in particular, also to the case of graphs. The uniqueness result is based on a 
very standard parabolic comparison principle. The general approach, based on the 
fact that the level sets of the distance functions have nonincreasing curvatures as 
the distance increases (as was exploited as early as in [39] in the viscosity setting), 
can quite probably be used in more general situations, and even maybe for motions 
which are not necessarily variational. However, it should need substantial adaption. 
For instance, if replacing the mobility m = <f>° in our approach by other (convex) 
functions is in principle easy (it is enough to consider, for the distance functions, 
the m-distance function instead of the (^“-distance), in the nonsmooth case it yields 
difficulties which still require further investigation. Indeed, if m is smooth and (f> 
is not, then it will not be true anymore that the level sets of the distance function 
have globally bounded curvature as the distance increases, so that Definition 12.11 
needs to be changed. It is not yet clear what assumption on (divz^ is then useful 
in order to be able to derive both existence and uniqueness. This is a subject for 
future study. 
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